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A key prediction of the 5*0(5) theory is the antiferromagnetic vortex state. Recent neutron scattering experiment 
on LSCO superconductors revealed enhanced antiferromagnetic order in the vortex state. Here we review theoretical 
progress since the original proposal and present a theory of static and dynamic antiferromanetic vortices in LSCO 
superconductors. It is shown that the antiferromagnetic region induced by the vortices can be greater than the 
coherence length, due to the light effective mass of the dynamic antiferromagnetic Ructuations at optimal doping, 
and close promixity to the antiferromagentic state in the underdoped regime. Systematic experiments are proposed 
to unambiguously determine that the field induced magnetic scattering originates from the vortices and not from the 
bulk. 
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SO{5) theory is an unified theory of antiferromag- 
netism (AF) and superconductivity (SC) in the cuprate 
superconductors |^,^. Although AF and SC phases seem 
to be detached in the experimental phase diagram, this 
theory predicts a direct transition between them as the 
chemical potential is varied. As a function of a physi- 
cal parameter, the direction of the SO{5) superspin vec- 
tor can be rotated smoothly. However, the chemical po- 
tential is hard to control experimentally. Doping in the 
AF/SC transition region is complicated by the chemical 
and structural inhomogeneities. Therefore, in the 5*0(5) 
theory, it was proposed that the interplay between these 
two phases can be investigated in a controlled way in the 
vortex state of the superconductor, where the external 
magnetic field can smoothly rotate the direction of the 
superspin vector. This theory makes the striking predic- 
tion that the core of the vortex is antiferromagnetic [Q , 
in sharp contrast to the metallic vortex core in a con- 
ventional superconductor. The AF vortex state (AFVS) 
can be viewed as a topologically non-trivial texture of 
the 5*0(5) superspin vector. Furthermore, it was explic- 
itly proposed that the AF vortex state can be detected in 
neutron scattering experiments as satellite peaks spaced 
by 27r/d, where d is the spacing between the vortex cores 
. These satellite peaks give a precise experimental def- 
inition of the AF vortex state, where static AF order 
coexists with the vortex state of the superconductor. 

A theory of the AF vortex core was developed by 
Arovas et al . This work introduced a method to study 
the vortex core based on a Schroedinger like equation for 
the magnetic excitations near the vortex core. Due to 
the 50(5) constraint on the AF/SC superspin vector 0, 
the vortex core acts like an attractive potential for the 
magnetic excitations. If a bound state exists, enhanced 
low energy magnetic fluctuation is localized near the vor- 
tex core. If the bound state energy goes to zero, static 
AF order is established inside the vortex core. Based on 
estimation of parameters, it was argued that the static 
vortex core is stable in the underdoped regime Q . These 
authors suggested that the AF vortex state can be de- 



tected in neutron scattering experiments, and predicted 
that the field induced antiferromagnetic moments should 
scale with the applied field, or the number of vortices in 
the system |3j. Soon after this work, Bruus et al Q| ana- 
lyzed experimental data in YBCO superconductors and 
argued that the vortex core in optimally doped materials 
should have enhanced dynamic AF order, with energy 
below the neutron resonance energy of AlmeV. Ogata 
Q performed extensive variational calculations on the 
t — J model and found direct evidence of the static AF 
vortex core up to doping level of (5 = 0.1. In a classical 
Monte Carlo calculation of the 5*0(5) model, Hu ^ ob- 
served the AF vortex state in the satellite peaks of the 
AF spin correlation function, which are spaced by 2it /d. 
Thermodynamic implications of the AF vortex state has 
been recently investigated by Juneau et al 0. Various 
groups argued that the AF vortex core is ex- 

pected to lead to strong suppression of the local density 
of states, consistent with the STM experiments Re- 
cently, Demler et al described the influence of the 
superflow kinetic energy and found that this led to an 
important contribution to the magnetic field dependence 
of the spontaneous moment in a phase with magnetic 
long-range order. Furthermore, these authors also de- 
veloped a theory for dynamic magnetism in the vortex 
lattice. This approach is reviewed in the same volume by 
Sachdev [||. 

On the experimental side, the search of AF vortex state 
showed promising results. Early on Vaknin et al 
found evidence of AF order in the vortex state of YBCO. 
Katano et al found enhanced AF fluctuations in the 
vortex state of LSCO. Most striking evidence of AF or- 
der in the vortex state is observed recently by Lake et 
al ||l6| , |l7| . In the vortex state of optimally doped LSCO 
superconductors, incommensurate magnetic fluctuations 
were found in the low energy region inside the spin gap, 
whose intensity scales with the number of vortices in the 
system. In the underdoped LSCO superconductors, en- 
hanced static AF order was discovered in the vortex state. 
In another related system, La2CuOi^y [|l8|,^, similar 
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field induced enhancement has been discovered. 

The purpose of this paper is to show that the theory of 
the AF vortex state can be applied straightforwardly 
to explain the recent neutron scattering measurements by 
Lake et al |16 1^. We discuss various energy and length 
scales in the problem, and show that experiments are 
performed in a regime where static and dynamic AF vor- 
tex state can be observed. In the original theory, only 
the case of the commensurate AF fluctuations were dis- 
cussed. As we shall see, the theory can be extended in a 
simple way to accommodate incommensurate AF fluctu- 
ations as well. Most importantly, we show that one can 
determine unambiguously whether the field enhanced AF 
fluctuation originates from the bulk or from the vortices, 
and show that present experiments already provides very 
strong evidence that the later case is realized. 

The starting point of our consideration is the effec- 
tive SO{5) theory of the magnetic fluctuations linearized 
in an inhomogeneous background of the SC condensate 
||,|. This theory can be cast into the following form 
of the effective Lagrangian £ = ^^i^^tnia)^ — H, where 
ma{x,t) is the local AF order parameter and 
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where uj — -^ApAg is the dispersion of the AF fluctua- 
tion in the absence of the external magnetic held. In the 
LSCO superconductors, we can approximate it with the 
following form: 
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Here is the spin gap energy. In optimally doped LSCO 
superconductors. As « 7meV. Qo is the wave vector of 
the incommensurate AF fluctuations in LSCO , which 
can be interpreted as dynamic stripes ||2l| , ^ . M* is the 
effective mass for the AF fluctuations, and we shall give 
an estimate of this parameter later. v{x) is the poten- 
tial due to the SC order parameter. In the SO (5) theory, 
there is a constraint between the AF and SC order param- 
eters in the form + = 1, where t/j is the SC order 
parameter. If we implement this constraint by a soft spin 
constraint, of the form —g{m'^ + IV'P) + w(to^ -I- 
we see that a repulsive coupling between the two order 
parameters is implied. This would remain true in the 
presence of a 5*0(5) symmetry breaking term wm'^\ip\'^, 
as long as w > —2u. The repulsion between the SC and 
the AF order parameters is the crucial physics for all field 
induced magnetic phenomena. If these two order param- 
eters were decoupled |2^], no such effects would exist. 
Linearizing with respect to this repulsive coupling in the 
bulk we obtain the spin gap A^. However, at the cen- 
ter of the vortex core, the SC order parameter vanishes, 
and the effective spin gap is lowered. Therefore, there is 
an effective attractive potenial v{x) for the AF magnetic 



fluctuations inside the vortex core. We assume that it 
takes the form 

v{x) ^^vq{x,Xi), vo{x,Xi)^-Ve (3) 

i 

where < Ag, ^ is the SC coherence length and Xi 
describes the center of the vortices. Equations (|^), (||) 
and determines the AF fluctuation spectrum in the 
vortex state. 

The description of the free propagation of the AF mag- 
netic fluctuation obtained from the linearized SO{5) the- 
ory is similar to that first obtained from analytically 
continuing a magnon triplet from a quantum disordered 
state, such as the spin Peierls state [^,^. However, the 
crucial point here is the attractive interaction between 
the AF magnetic fluctuations with the vortices, which is 
flrst obtained from the SO{5) theory. In the original the- 
ory, a commensurate form of the dispersion is assumed. 
However, it was later found that competing interactions 
in the S0{5) theory could produce dispersion in the form 
of (|^), with a general incommensurate wave vector Qo 
p6[ . For our current discussion, it is easy to see that Qo 
can be trivially gauged away from the problem. If we 
perform a phase transformation 
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we see easily that the quadratic Hamiltonian can be diag- 
onalized by the following Schroedinger like equation for 
na{x): 
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with energy given by As -I- A. This equation was flrst 
introduced in equation (10) of ref. and forms the basis 
for subsequent studies on AF order in the vortex state. 

Equation (|^) describes the motion of an effective quan- 
tum particle with mass M* /2 in a periodic array of at- 
tractive vortex centers. The important physics here is 
that in D=2, at least one bound state forms inside the 
attractive potential. This bound state for the collective 
AF fluctuation is not to be confused with the bound state 
of quasi-particles inside the vortex cores. Let us first con- 
sider the simple case where only one bound state forms 
inside the attractive potential, with energy Eb- The lo- 
calization length of the bound state can be estimated to 
be 
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In the limit where I < d, the band width of the vortex 
band can be estimated to be 



t - Ve^ 



(7) 
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The three length scales, the SC coherence ^, the AF mag- 
netic localization length the inter- vortex spacing d, and 
the four energy scales, the spin gap Ag, the attractive 
potential V, the bound state energy Eb and the vortex 
band width t are illustrated in Fig. 1. The concept of 
a static or dynamic AF vortex core is a sharply defined 
concept. If there is no bound state, then there is no 
sense in which one can define the concept of AF vortex 
core. However, if bound states exists inside the spin gap, 
the vortex core develops dynamic AF ordering, with well 
defined ordering energy scale As — Eb and well defined 
ordering length I. One can view I as the size of the vortex 
and due to the finite size quantization, the AF moment 
fluctuates dynamically, with time scale h/(As — Eb)- In 
the regime of weak field, I << d, the bandwidth of the 
AF vortex band can basically be neglected. With in- 
creasing field, the bound state energy Eb decreases, and 
the inter-vortex coherence t increases as well. When the 
bottom of the vortex band energy touches zero, the static 
AF vortex state is obtained. 

With this basic preparation we can now discuss the 
experimental situation in the optimally doped LSCO su- 
perconductors, and present the central argument of this 
paper that the field induced scattering originates from 
AF fluctuations localized around the vortex cores. Ex- 
periments are done in the weak magnetic field regime 
where the volume fraction associated with the vortices is 
extremely small, / ^ {(,/d)'^ < 10%. On the other hand, 
the field induced magnetic scattering is centered around 
AmeV, well below the spin gap energy of = 7meV. 
If the field induced AF magnetic fluctuations were ex- 
tended, scattering by the vortices can only cause an en- 
ergy shift proportional to V{^/d)'^, which vanishes in the 
limit of low vortex volume density. See Fig. ^ In the 
experimental regime oi f ^ 10%, the energy shift of ex- 
tended states would be of the order of AE — Vf < 
Agf < O.lmeV, far less than the experimentally observed 
value of AE ~ 3meV. Therefore, the Lake et al experi- 
ment |l^ already shows that it is extremely unlikely that 
the field induced scattering originates from the bulk. On 
the other hand, if the field induced AF magnetic fluctu- 
ations are localized around the vortices, the energy shift 
is finite in the limit / — > 0. See Fig. ^. If we assume 
one bound state per vortex, the field induced intensity 
/ is proportional to /. Therefore, this picture predicts 
that the field induced intensity is directly proportional to 
H/ Hc2 ■ These are the precise signatures of the AF vortex 
core in the limit of low vortex volume density. In Fig. 2 
of ref. 1 16 1, one sees clearly that the continuum above the 
spin gap is little changed by the field, but new spectral 
weight is introduced with a centroid well detached from 
the continuum. The experimental fact that the centroid 
of the field induced signal is well below the spin gap and 
the intensity is proportional to / provides strong evidence 
that the field induced scattering is localized around the 
vortices. A more systematic analysis should plot both 



the energy centroid and the intensity of the field induced 
signal as a function of /. Since t depends exponentially 
on d, it vanishes faster than any other energy scales in 
the limit of / — > 0. Therefore, if Eb — > finite and I (x f 
in the extrapolated limit o/ / — > 0, unambiguous evidence 
in support of the AF vortex core can be established. The 
/ ^ limit is an useful tool to label the bands, since once 
/ becomes finite, there is no rigorous distinction between 
the extended and localized bands. 

Considerations of the energy of the field induced scat- 
tering leads to the conclusion that the they must be lo- 
calized near the vortex cores. Let us now fix the energy 
at As — Eb and consider the spatial distribution of the 
AF fluctuations. The three length scales can be easily es- 
timated from the experiments, giving ^ 4ao, I ~ 20ao 
and d ~ 40ao. Here I is estimated from the width of 
the momentum distribution function at flxed energy of 
AmeV. Theory predicts ||l|,^,^ that the momentum dis- 
tribution function for the n{x) fleld consists of a series of 
satellite peaks spaced by G = 2iT/d, centered around the 
main peak at g 0, where the intensity of each satel- 
lite peak is covered by an envelope function of the width 
l/l. See figure (^. Neutron scattering couple to the 
m(x) field, by applying the transformation (^) we obtain 
the same pattern, but with the center of the envelope 
function shifted to incommensurate wave vector Qo- See 
figure (^. Since the momentum shift, or gauge, transfor- 
mation applies equally to both extended and localized 
states, the AF fluctuation localized near the vortex core 
is associated with the exactly the same incommensurate 
wave vector Qq as the AF fluctuations in the bulk. 

Now we need a quantitative argument on why the ob- 
served AF coherence length I 20ao seems longer com- 
pared to ^. From equation (^, we see that this could 
be due to the light effective mass AI* . But M* could be 
independently measured in the neutron scattering exper- 
iment in the absence of the applied field, by fitting the ex- 
perimentally measured AF dispersion with tUp = y^ApAg 
given in equation . This provides an independent con- 
sistency check on our theory. A small effective mass 
translates into steep dispersion in Up versus p. Exper- 
imentally, the incommensurate peak of the momentum 
distribution function for various energies does seem to 
shift, implying a very small effective mass. If we as- 
sume that the incommensurate peak does not shift over 
an energy range of 8meV, and momentum shift less than 
(lOflo)"^ can not be detected, we obtain a value for M* 
which is roughly consistent with the coherence length 
I ^ 20ao. It would be highly desirable to improve the ex- 
perimental accuracy and quantitatively check the relation 
between the inelastic neutron scattering in the absence 
of the field, bound state energy and the coherence length 
in the vortex state. 

As the external magnetic field increases, d decreases 
and the intervortex coherence t increases. For d I, 
it ceases to take the form t ~ Ve~'^^^. A similar effect 
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arises when the system is more underdoped, where both 
^ and I increases, resuhing in the increase of t. The local 
bound state near each vortex now form a band, with band 
width t. As a function of the applied field, the center of 
the band shifts downwards, and the bandwidth increases. 
Within Born scattering, the shift of the band minimum 
can be estimated to be 



Emin = Eb - aVf 
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where a is a dimensionless constant. When the bottom of 
this vortex band touches zero, the system makes a phase 
transition into a AF vortex state as predicted in 
However, Demler et al |12 observed that in this regime, 
the superflow makes an important correction to the field 
dependence, and find a logarithmic correction As 
argued in previous paragraphs, the static AF moments 
are not localized in a region of the size of ^, but of the 
size of I, which is now comparable to d. It is important to 
point out here that due to the repulsion between AF and 
SC order parameters, the static AF order near the vor- 
tices is always larger than the AF order in the bulk. How- 
ever, these quantitative corrections do not change the 
basic periodic structure. As pointed out in ref. ||i|,|6|,p^ , 
this periodic pattern can be detected directly in neutron 
scattering experiments as satellite peaks. One should ob- 
serve the momentum distribution functions displayed in 
Fig. 2 and 3, for quasi-elastic scatterings. These satellite 
peaks provide direct experimental evidence for the static 
AF order in the vortex state. Recent neutron scatter- 
ing experiments in underdoped LSCO superconductors 
revealed long ranged static incommensurate AF order in 
the presence of an external magnetic field where 
the field induced signal is proportional to H/Hc2, or the 
number of vortices in the sample. It is highly desirable to 
refine the momentum resolution in the experiments and 
detect the satellite peaks directly. 

Now we are in a position to understand why relatively 
moderate field Hafvs can induce a transition into a 
statically ordered AF vortex state. If there were no lo- 
calized states near the vortices in the / — *■ limit, the 
transition to a statically ordered state arises from the 
band of extended states, whose bottom touches zero. 
However, if the band of localized state exists below 
the extended band, its bottom will always touch zero 
first. Since bands are not expected to cross as the 
external magnetic field is varied, we can roughly esti- 
mate the ratio between these two critical fields to be 
Hafvs (bound)/ Hafvs {extended) ^ i^s — Eb)/^s ~ 
1/2 < 1. In other words, transition due to the localized 
band always preempts the transition due to the extended 
band. From this argument, we reach the important con- 
clusion that the relatively moderate field can induce a 
transition into a statically ordered AF vortex state pre- 
cisely because the AF order originates from the vortices. 

In above discussions we mostly focused on neutron 
scattering experiments in LSCO. However, the theory is 



generally applicable to all high Tc materials. In recent 
NMR experiments , enhanced AF fluctuations around 
the vortex cores have been detected in the YBCO system. 
In a recent STM experiment, Hoffman et al |^ detected 
CDW-like order around the vortex core in the BSCO sys- 
tem. In both cases, the magnetic field induced order are 
well localized around the vortex core, with a localization 
length greater than the vortex core itself. The important 
point to be stressed here is that the localization length 
is finite, in accordance to the general discussions out- 
lined above. It would be highly desirable to extend the 
neutron scattering experiments under the magnetic field 
already performed by Dai et al |2^ on YBCO systems to 
higher field, and to look for the dynamic AF bound state 
around the vortex core, with bound state energy of few 
meV below the neutron resonance energy. 

Finally, we would also like to comment on a new order 
parameter which emerges in the AF vortex state. In ref. 

one of us introduced a general constraint between 
the SO{5) superspin order parameters ria and the SO{5) 
symmetry generators Lab, which takes the form: 
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In regions of space where both the AF order parame- 
ter n2,3,4 and the SC order parameters ni^s coexist, the 
SO{5) constraint (||) automatically implies that the cor- 
responding TT operators ^12,13,14 and ^25,35,45 would ac- 
quire non-zero expectation values as well. Therefore, in 
the region ^ < r < I around the vortex core, where both 
AF and SC orders coexist, the order associated with the 
TT operators in the SO{5) theory exists as well. In a sub- 
sequent paper, we shall fully explore the experimental 
consequence of this observation. 

In conclusion we have reviewed theoretical progress 
since the original proposal of AF order in the SC vor- 
tices [0J|], and presented a semi-quantitative theory of 
the static and dynamic AF order in the vortex state of 
LSCO superconductors. It is argued that the field in- 
duced signals found in experiments arise from the dyan- 
mic and static AF order in the superconducting vortices. 
We showed that the intermediate length scale I describes 
the localization of AF order near SC vortices. This length 
scale is large in the LSCO materials because the effective 
mass for the incommensurate AF fluctuations in the ab- 
sence of the field is anomalously small. However, this 
length scale is finite in the limit of low vortex density. 
Whether the field induced scattering originates from the 
bulk or from the vortices can be systematically distin- 
guished in neutron scattering experiments, where the en- 
ergy shift of the field induced scattering is plotted against 
H, in the limit of low H . If the energy shift remains fi- 
nite in the zero field limit, an unambiguous case for the 
AF vortices can be established. Satellite peaks are pre- 
dicted in the momentum distribution functions for both 
quasi-elastic and inelastic scatterings, revealing the static 
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and dynamic AF order of vortices. An exact argument 
based on gauge transformation is presented to show that 
the AF order locahzed around the vortices has the same 
incommensurate wave vector as the bulk AF fluctuations. 

The theoretical prediction and experimental discovery 
of the AF vortex state offers crucial insights into the mi- 
croscopic mechanism of high Tc superconductivity. Vari- 
ous theoretical proposals differ only in the nature of the 
intermediate state between the AF and the SC state. If 
doping or chemical potential were the only route connect- 
ing these two phases, the problem of chemical inhomo- 
geneities will always obscure the underlying physics. The 
external magnetic field provides a clean alternate route 
to connect these two phases and a new route to attack 
this problem. Unlike the case of chemical doping, exper- 
iments, coupled with the scaling arguments presented in 
this work, can precisely determine the real space loca- 
tions of the AF signals. While the 5*0(5) theory predicts 
an AF vortex state 0, other theories predict hc/e vortex 
pO| , visons 1^ and staggered flux phase |Q. Therefore, 
the nature of the vortex state is important in distinguish- 
ing among the various theoretical proposals. 
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FIG. 1. Illustration of the energy and length scales in 
the problem. As is the spin gap energy, V is the attractive 
potential due to SC vortices, Eb is the bound state energy 
and t is the bandwidth of the AF vortex band. ^ is the SC 
coherence length, I the localization length of AF fluctuations 
around the vortices and d is the inter-vortex spacing. 
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FIG. 2. If no bound states exist, energy of extended states 
can only shift by an amount of the order oiV{£_/d)'^ , which is 
vanishingly small in the limit of d ^ 0. 
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FIG. 3. In the presence of bound states, the energy shift 
is finite, while the intensity due to the bound state band is 
proportional to {^/d)^ ^ H/Hc2- 




Qo 2G Qo G Qo Q0 + G0Q0 + 2G q 
FIG. 4. At energy E — As ~ Eb for optimal doping and 
E = Q for underdoping, the momentum distribution func- 
tion consists of a series a delta function peaks, spaced by 
G = 2iT/d. The envelope of these peaks forms a broader peak 
centered around the incommensurate wave vector Qo, with 
width l/l. 




FIG. 5. A more explicit illustration of the satellite peaks 
near the four incommensurate spots. The inset shows the 
vortex lattice in real space. 
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